Creation of pure multi-mode entangled states in a ring cavity 
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Practical schemes for creation of multi-mode squeezed (entangled) states of atomic ensembles located inside a 
high-Q ring cavity are discussed. It is assumed that the cavity is composed of two degenerate mutually counter- 
propagating modes that can simultaneously couple to the atomic ensembles with the same coupling strengths. 
The ensembles are composed of ultra-cold atoms which are modeled as four-level systems driven by two laser 
fields, both co-propagating with one of the cavity directions. We illustrate a procedure that constructs multi- 
mode squeezed states from the vacuum by a unitary transformation associated with the collective dynamics of 
the atomic ensembles subjected to driving lasers of a suitably adjusted amplitudes and phases. The lasers pulses 
together with the cavity dissipation prepare the collective modes in a desired stationary squeezed state. 
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I. INTRODUCTION 

The problems of generation and practical applications of 
squeezed states of the electromagnetic field have been under a 
continuous interest since the early days of quantum optics Q|- 
yfl. Squeezed states, predicted in many phenomena of quan- 
tum optics and laser physics, are an example of non-classical 
states of the field which are characterized by a reduction of 
the quantum fluctuations in one of the field quadrature com- 
ponents below the usual vacuum level. The best description of 
how many different processes can exhibit squeezing was given 
by Krzysztof Wodkiewicz, who comments "it is very easy and 
rather trivial to give many examples of quantum-mechanical 
states that lead to squeezed quantum fluctuations" 0]. The 
ability to achieve a good quality squeezed light with signifi- 
cantly reduced fluctuations has not only resulted in the devel- 
opment of new areas of research but also provided a possibil- 
ity of testing many fundamental ideas of quantum physics Jj- 
Ht|. Apart from the development of new theories of the in- 
teraction of quantum systems with the electromagnetic field, 
a number of practical applications of squeezed light have 
been proposed and realized ranging from optical communica- 
tions Hi, gravitational- wave detection ifToll to newly emerging 
quantum technologies such as quantum information process- 
ing fTD - Tnll. quantum teleportation ITii [T3|. quantum cryp- 
to graph v Ilfj|,[l7tl and quantum computing Il8lfl9ll . 

Recent investigations of the continuous-variable quantum 
information have also generated a great deal of interest in 
squeezed states of light i20ll . A particular attention has been 
paid to generation of multipartite continuous variable entan- 
gled states of atomic ensembles form multi-mode squeezed 
light produced, for example, in a non-linear process of para- 
metric down conversion l2lU24ll . The reason for using atomic 
ensembles is twofold. On the one hand, atomic ensembles are 
macroscopic systems that are easily created in the laboratory. 
On the other hand, the collective behavior of the atoms en- 
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ables to achieve almost a perfect coupling of the ensembles to 
external squeezed fields without the need to achieve a strong 
field-single atom coupling. 

Particularly interesting are schemes based on a quantum 
reservoir engineering, where squeezed states are generated in- 
side the atomic ensembles by internal dynamical processes 
rather than being injected into the ensembles from external 
sources. The squeezed states are generated by a suitable driv- 
ing that creates an effective multi-wave mixing processes to 
occur in the interaction of the fields with the atoms. For ex- 
ample, Rice and Carmichael l25ll have shown that the fluores- 
cence spectrum of single two-level atoms driven by a weak 
coherent laser field can be narrowed below the quantum limit 
giving a subnatural linewidth of the emitted field. The nar- 
rowing has been explained as resulting from squeezing in the 
fluctuations of the atomic dipole moments induced by the in- 
teraction of the atoms with the laser field. Narrowing of the 
spectral line below the quantum limit has also been predicted 
in the output field of a high-Q microcavity engineered within 
a photonic crystal and containing two-level atoms driven by 
a strong laser field l26ll . In this case, the spectral narrowing 
also results from squeezing, but the mechanism is different. 
The squeezing is produced by a non-linear mixing process re- 
sulting from a non-linear interaction of the cavity field with 
the dressed atom. 

In a recent paper, Parkins et al. ll27ll have proposed a 
scheme for preparing single and two-mode squeezed states in 
atomic ensembles located inside a high-Q ring cavity. The 
scheme, which is a generalization of the Guzman et al. ll28ll 
work to four-level atoms, is based on a suitable driving of 
the atomic ensembles with two external laser fields and cou- 
pling to a damped cavity mode that prepares the atoms in a 
pure squeezed (entangled) state. Similar schemes have been 
proposed to realize an effective Dicke model operating in the 
phase transition regime l29ll . to create a stationary subradiant 
state in an ultracold atomic gas lT30ll . This approach has re- 
cently been proposed as a practical scheme to prepare trapped 
and cooled ions in pure entangled vibrational states [31j| and to 
prepare four ensembles of hot atoms in pure entangled cluster 
states JH. 

The scheme proposed by Parkins et al. l27ll assumes that 
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atomic ensembles are effectively coupled to a single or two 
modes co-propagating with the driving fields. In fact, a prac- 
tical ring cavity is composed of two degenerate mutually 
counter-propagating modes that can simultaneously couple to 
the atomic ensembles with the same coupling strengths ll33l - 
l35ll . With the counter propagating mode included the situ- 
ation is somewhat more complicated, since the presence of 
the counter-propagating mode may result in a mode coupling 
which depends on atomic positions ll36ll . 

In this paper, we propose to generalize the Parkins et al. 
model to the case of a high-Q ring cavity with two mutu- 
ally counter-propagating modes, as illustrated in Fig. Q] In 
this two-mode configuration, the photon redistribution can 
occur between the counter-propagating modes which, as we 
shall see, will result in multi-mode squeezing. We propose a 
simple procedure to create multi-mode squeezed (entangled) 
states in the practical scheme involving cold atomic ensembles 
composed of four-level systems driven by sequences of laser 
pulses of specifically chosen Rabi frequencies and phases. 
The procedure constructs squeezed states from the vacuum by 
a unitary transformation associated with the dynamical pro- 
cess determined by the master equation of the system. With 
this more practical model, we illustrate how to create single 
and two-mode squeezed states in a single ensemble, and a 
four-mode squeezed state between two atomic ensembles. 




FIG. 1 : A schematic diagram of a two-mode ring cavity containing 
two ensembles of cold atoms trapped along the cavity axis. Both 
cavity modes are damped with the same rate k. The driving laser 
fields are injected through the cavity mirrors and co-propagate with 
one of the cavity modes. 



a high-Q ring cavity. Each ensemble is composed of TV 
identical multi-level atoms located at different positions Xj n 
and modeled as four-level systems with two ground states 
\0j n ), |ljn) and two excited states \uj n ), \sj n ), where the sub- 
script jn labels jth atom of the nth ensemble, see Fig. [2] 
In practice such a four-level system could correspond to an 
F = 1 -<-> F' = 1 transition as occurs in 8 'Rb atoms. The 
cavity is composed of three mirrors that create two mutu- 
ally counter-propagating modes, called clockwise and anti- 
clockwise modes, to which the atoms are equally coupled. 
The cavity modes are degenerate in frequency, i.e. uj + = 
uj- = oj c . Each ensemble is driven by two laser fields in- 
jected through the cavity mirrors and co-propagating with one 
of the cavity directions. The frequencies ljl u and ljl s of the 
laser fields are matched close to the cavity frequency, so that 
the wave numbers of the laser fields k u and k s are approxi- 
mated by k u w k s k, k. We assume that the laser fields and 
the cavity frequencies are far detuned from the atomic transi- 
tion frequencies. This will allow us to adiabatically eliminate 
the upper states of the atoms and obtain an effective dispersive 
type interaction of the driving lasers and the cavity modes with 
the atomic transitions. 
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FIG. 2: Outline of the four-level system driven by two highly de- 
tuned laser fields and simultaneously coupled to cavity modes. The 
laser fields of the Rabi frequencies fl un and Q sn drive the atomic 
transitions |lj n ) — > \ u jn) and \0j n ) — > \sj n ), respectively. The 
cavity modes are coupled to atomic transitions \lj n ) — > \sj n ) and 
\0j n ) — s- \ujn) with the coupling strengths g sn and g un , respec- 
tively. 



Construction of pure squeezed states via unitary transfor- 
mations, which we consider here, is a fascinating research 
topic in which Krzysztof Wodkiewicz was interested for many 
years. In a series of publications, he has explored similarities 
between the transformations and simple Lie algebras, and also 
investigated connections between squeezing and other quan- 
tum effects, such as the EPR paradox, photon antibunching 
and entanglement l37l - l43ll . 



II. EFFECTIVE HAMILTONIAN 

The physical system considered consists of two one- 
dimensional atomic ensembles trapped along the axis of 



The Hamiltonian of the system, in the rotating-wave ap- 
proximation, has the following form 



H = H + H AL + H AC , 



(1) 



where 
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H = fko c (a\a + + aLa_) + h^^{u> u \uj n )(uj n \ 

71=1 j = l 

+ uj s \s jn )(s jn \ +cJi|lj„)(lj- n |} (2) 
is the free Hamiltonian of the atomic ensembles and the cavity 
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modes, 

Hal = ^^X!{^n(^-„)e-^» + ^^ t |u in }(l in | 

n=l j=l 

+ O s „(^ n )e-^^ s+ ^ t |s j „)(0 i „| +H.c.} (3) 



2 AT 



is the interaction Hamiltonian between the atoms and the driv- 
ing fields, and 
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H A c = h J2J2i [9tn( x jn)a+ + g un {x jn )aJ\ K-„)(0,- n | 

n=l 3 = 1 

+ [9tn( x in)a+ + g~ n (x jn )a-] \s jn ){l jn \ + H.c.} (4) 



is the interaction Hamiltonian between the atoms and the two 
cavity modes. 

Here, a± and a± are the annihilation and creation oper- 
ators associated with the two counter-propagating modes of 
the cavity; clockwise (+) and anti-clockwise (— ) propagat- 
ing modes, and k is the wave number of the cavity modes. We 
have denoted the energies of the atomic levels by huii (i = 
1, u, s) and have set the energy of the ground state \0j n ) equal 
to zero. The parameters Vl un (xj n ) and il sn (xj n ) are the po- 
sition dependent Rabi frequencies of the driving laser fields, 
and (f> un , 4> S n are their phases. The coupling constants of the 
atomic transitions to the cavity fields, gt n ( x jn) and gf n (xj n ), 
are also dependent on the atomic position. In what follows, we 
will use the plane traveling wave representation for the laser 
fields and the cavity modes, in which 



9un( x jn) — gunC 



9sn{ X jn) — 9snC J ™ i (5) 



where we have assumed that the coupling constants g un 
and g sn are the same for the two cavity modes. This is accept- 
able if the modes have the same polarization and geometry, 
which is feasible with the current experiments 03143511 . 

We now make an unitary transformation and few standard 
approximations on the Hamiltonian (Q]i to eliminate the ex- 
plicit time dependence and the excited states to neglect atomic 
spontaneous emission and to obtain an effective two-level 
Raman-coupled Hamiltonian. Then, we will make the trans- 
formation of the atomic operators into the field (bosonic) rep- 
resentation. These operations will be done in the following 
steps. In the first step, we make the unitary transformation 
U = cyLY>{iH' t/K), with 
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+ H^2^{(UJ U + U) 1 )\u jn )(Uj n \ 
n=l j=l 

+ULs\Sjn)(Sjn\ + ^1 1 Ijn) (Ij'n 1} , 



and assume that the laser frequencies satisfy the resonance 
condition Wi s — W£.„ = 2uii. This specifically chosen trans- 
formation allows us to eliminate the explicit time dependence 
of the Hamiltonian. 



In the second step, we introduce detunings of the laser fields 
from the atomic transition frequencies 



A u = u> u - (u> Lu A s = u a - u Ls , 



(7) 



and assume that the detunings are much larger than the Rabi 
frequencies and the atomic spontaneous emission rates. This 
allows us to perform the standard adiabatic elimination of the 
atomic excited states and obtain an effective two-level Hamil- 
tonian of the form 
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where <5 C = uj c — (w£ S — u>i) is the detuning of the cavity 
frequency from the Raman coupling resonance, 
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imkxj 



are position dependent collective atomic operators, and 
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are the coupling strengths of the effective two-level system 
to the cavity modes. The three collective atomic operators, 

Jak i^2k anc ^ J -2k ar ^ se naturally for the position depen- 
dent atomic transition operators and appear in a cavity with 
two mutually counter-propagating modes. In the case of a 

(n) 

single-mode cavity, the Hamiltonian involves only the J^ k 
operator (22I HI 113] ■ The operators crj' l) , oj n)t and a ( z f are 
the standard Pauli spin— 1/2 operators, which satisfy the well 
known commutation relations [aj n \a < f 1 ^] = 2aJ^SjiS mn . 

However, the collective operators J^lt ^mk anc ^ c '° not 
in general satisfy the angular momentum commutation re- 
lations. The reason is in the presence of the phase factors 
exp(imkxj n ), so that the commutation relations are satis- 
fied only in the small sample limit of kxj n < 1, at which 
(6) cxTp(imkxj n ) ~ 1, 

To avoid unessential complexity, we choose the frequencies 
of the driving lasers and the cavity field such that 
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and find that after this simplification, the Hamiltonian <(8j re- 
duces to 



H = 



T J 2k a - 



+ fee"*" (j&M + J^V ) + H.c.} . (12) 

This equation is in the form of a non-RWA Hamiltonian of 
two extended atomic ensembles independently coupled to two 
counter-propagating cavity modes. The parameters of the 
Hamiltonian are a function of the detunings and Rabi fre- 
quencies of the two highly detuned laser fields co-propagating 
with the clockwise cavity mode and thus could be controlled 
through the laser frequencies and intensities. In the small 
sample case of kx ]n <C 1 and under the single-mode ap- 
proximation, the Hamiltonian simplifies to the standard Dicke 
model 00. 

In the final step we shall reformulate the Hamilto- 
nian (fT2l in terms of bosonic variables by adopting the 
Holstein-Primakoff representation of angular momentum op- 
erators 1 46]. In this representation, the collective atomic oper- 
ators, , J^l and Jg are expressed in terms of annihila- 



tion and creation operators and 1 of a single bosonic 
mode. Provided the atoms in each ensemble are initially pre- 
pared in their ground states {|0 jn }}, and taking into account 
that due to large detunings of the driving fields, the excitation 
probability of each atom is low during the laser-atom-cavity 
coupling, i.e., (c^™ ) ~ — 1/2, the collective atomic operators 
can be well approximated by 



Ji"2 = Vnc. 



(n) 
ik ' 



Jin) _N 
z ~ 2 ' 



(13) 



where 
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(n) 
mk 
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c (n) e<tn te jB m = ,±2, (14) 



(™) 



are collective bosonic operators with the operators b 
and b ™ obeying the standard bosonic commutation relation 

rt(") .(m)ti s s 

[Oj ,o t \ - 0j£0 nm . 

Note that the collective bosonic operators do not in general 
commute, i.e. 



fl(n) y-f(f»')t 

mk ' m'k 



— cxp[i(m 
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m )kxj]6 n 



(15) 



Hence, the modes are not independent of each other. 

However, suppose that the atoms are in a chain with a con- 
stant separation d between the adjacent atoms, which is much 
smaller than the cavity wavelength, i.e., kd <C 1. Then the 
position of the jth atom can be expressed as Xj = (j — l)d, 
and the commutation relation (TT~5T > simplifies to 



£<(«) /~<( n ')t 
mk ' m'k 



exp[i(m — m')kL] — 1 
i(m — m')kL 



(16) 



where L — Nd is the length of the atomic ensemble. 

If the length of the atomic medium is much larger than the 
cavity wavelength, L ^> A, we obtain 
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which shows that in the limit of small separations between the 
atoms and long atomic chains, the the collective bosonic oper- 
ators are orthogonal to each other. One can argue that at small 
separations between the atoms, assumed in the above deriva- 
tion, the direct dipole-dipole interaction between the atoms is 
large and should be included into the calculations. However, 
we assume that the atoms are initially prepared in their ground 
state and owing the fact that the laser and the cavity mode fre- 
quencies are far from the resonance between the ground and 
excited states of the atoms, no transition dipole moments are 
induced between the atomic states. Thus, the direct dipole- 
dipole interaction between the atoms can be ignored. 

(n) (n) 

Using the collective bosonic operators C ok and Cjj_ 2fe we 
can rewrite the effective Hamiltonian ( fT2l in the following 
form 



n=l 

+ /Sune-*- ( + nC^l) ol + H.c.} .(18) 



£H e -i(0 s 



(19) 



The important property of the bosonic representation is the 
fact that effective Hamiltonian of ensembles of cold atoms 
trapped inside a ring cavity can be expressed now as the in- 
teraction between the cavity modes and three orthogonal field 

(n) 

modes; a collective mode C^ k solely coupled to the cavity 
mode o+, which co-propagates with the driving lasers, and 

(n) 

two modes C±2 k that are solely coupled to the cavity counter- 
propagating mode a_. 

The Hamiltonian ( TT~8T > holds for the laser fields co- 
propagating with the clockwise mode only. Following the 
same procedure as above, we can easily show that in the case 
of the driving fields co-propagating with the anti-clockwise 
mode a_, the effective Hamiltonian takes the form 

2 

H e = £ {(^CW +/3 fln e- i *"C& )t ) al 

n=l 

+ (pune-^C^l+psne-^cif^ 4 + H.c.} (20) 

We see the complete symmetry between the two cases that re- 
versing the direction of the propagation of the laser fields from 
clockwise to anti-clockwise is equivalent to the exchange of 
a + <-> a_ and k — > — k in the Hamiltonian ( TT8l . 

Our objective is to prepare the atomic ensembles in a de- 
sired pure squeezed vacuum state. To achieve it, we con- 
sider the evolution of the system under the effective Hamil- 
tonian ( fT8l including also a possible loss of photons due to 
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the damping of the cavity mode. This is the only damp- 
ing which we will consider as we have already eliminated 
or minimized spontaneous emission by choosing large detun- 
ings of the driving lasers. The spontaneous emission rate 
due to off-resonant excitation is estimated at the rate j c g = 
j(7/27r)(Q u ( s )/A^/ s s that with typical experimental values 
of 7 = 6 MHz for a rubidium atom and Sl u ( s )/A ll ( s ) = 0.005 
gives 7 c ff ps 40 Hz. The estimated value for j e g is signif- 
icantly smaller than k predicted for a cavity of the finesse 
F = 1.7 x 10 5 . 

With the cavity damping included, the properties of the sys- 
tem are determined by the density operator p whose the time 
evolution is governed by the master equation 



p = -i[H e ,p] +C c p, 



(21) 



C c p = -k haipal - a\a,ip - pa\a^j , (22) 



where 



i=± 

is an operator representing the damping of the cavity field 
modes with the rate k. 

In the following sections we will demonstrate how to gen- 
erate on demand multi-mode squeezed vacuum states in en- 
sembles of cold atoms located inside a two-mode ring cavity. 
We propose a simple procedure which constructs squeezed 
states from the vacuum by a unitary transformation associated 
with the realistic dynamical process determined by the master 
equation (|2TT >. In the procedure the atomic ensembles, ini- 
tially in the ground state, are subjected to different sequences 
of laser pulses of specifically chosen Rabi frequencies Vl un 
and fl sn and phases <fi un and <f> sn . The laser pulses together 
with the cavity dissipation, that occurs with the rate k, prepare 
the ensembles in a desired stationary state. We will assume in 
all our considerations that the initial ground state of the atomic 
ensembles corresponds to all the atoms being in their ground 
states |0in). 



III. CREATION OF ONE AND TWO-MODE SQUEEZED 
STATES IN A SINGLE ATOMIC ENSEMBLE 

First we illustrate a procedure which constructs single and 
two-mode squeezed states in a single (n = 1) ensemble of 
cold atoms located inside a two-mode ring cavity. The proce- 
dure constructs the squeezed states simply by acting on the 
vacuum with unitary operators called single and two-mode 
squeezed operators, respectively. The squeezed operators are 
defined as |47il 



So(£o) 



exp 



S±*(6) - exp (£C 2fc C_ 2fc - 



-2k j ' 



(23) 



where £o an d £i are complex one and two-mode squeezing pa- 
rameters, respectively. The squeezed operators can be easily 
associated with evolution operators for the effective Hamil- 
tonians ( fT8l and ( f20l i. which are already expressed in terms 



of the collective bosonic operators and C^k- 



The 



squeezing parameters are then given in terms of the coupling 
strengths (3 un and f3 sn and the phases 4> un , 4> S n, so they can 
be adjusted and controlled by the driving laser fields. 

The construction procedure is done in two steps. In the first 
step, we adjust the driving lasers to propagate in the clockwise 
direction, along the cavity mode a+. In this case, the dynam- 
ics of the system are determined by the Hamiltonian ( fT8l ). We 
then send series of laser pulses of phases u i = <f) s \ = 
and arbitrary Rabi frequencies fl u i and f2 s i, but such that 
Aa > Psi- With this choice of the parameters of the driv- 
ing lasers and under the unitary squeezing transformation 



with 



So(-£o)S±k(-£i)pSo(£ )S±k(£i) = p, 



the master equation (f2TJ becomes 



— p= -i[H e ,p] +£ c p, 



(24) 



(25) 



(26) 



where 



H e = So(-£o)S±fc(-£l)ffeSote>)S±*(6) 



H.c 



(27) 



It is seen that under the squeezing transformation, the Hamil- 
tonian represents a simple system of two independent linear 
mixers, where the collective bosonic modes Cofc and C-ik lin- 
early couple to the cavity modes a + and a_, respectively. The 
mode C-2k is decoupled from the cavity modes and therefore 
does not evolve. In other words, the state of the mode C_2fc 
cannot be determined by the evolution operator for the Hamil- 
tonian ( |27] >. The important property of the transformed system 
is that the master equation ( f26b is fully soluble, i.e. all eigen- 
vectors and eigenvalues can be obtained exactly. Hence, we 
can monitor the evolution of the bosonic modes towards their 
steady-state values. Since we are interested in the steady-state 
of the system, we confine our attention only to the eigenvalues 
of Eq. ( 126*1 ). which are of the form 



V± 



(28) 



Evidently, both eigenvalues have negative real parts which 
means that the system subjected to a series of laser pulses up 
to a short time t will then evolve (decay) to a stationary state 
that is a vacuum state. Thus, as a result of the interaction given 
by the Hamiltonian (f27b . and after a sufficiently long evolution 
time, the modes a±, Cofc and C2k will be found in the vacuum 
state, whereas the mode C_2fc will remain in an undetermined 
state. The state of the mode C_2fc will be determined in the 
next, second step of the preparation process. 

In order to estimate the time scale for the system to 
reach the steady-state, we see from Eq. ( f28l ) that as long as 
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\J fiti Pli > the time scale for the system to reach 
the steady state is of order of ~ 2/k. Thus, as a result of the 
cavity damping the system, after a sufficient long time, will 
definitely be found in the stationary state. 

In summary of the first step of the preparation, we find that 
in the steady-state, the density matrix representing the state of 
the transformed system is in the factorized form 



p(r ~ 2/k) = p v <g) p c _ 



(29) 



where 



Pv = |0 Q+ ,0 a _,0c Ofc ,0c 2 J(0 Q+ ,0 Q _,0 Cofc ,0c 2 J (30) 

is the density matrix of the four modes prepared in their vac- 
uum states, and pc _ 2k is the density matrix of the mode C-2k 
whose the state has not been determined in the first step of the 
procedure. The ket |0 a+ , a _ , 0c O k > 0c 2fc ) represents the state 
with zero photons in each of the modes. 

Thus, we are left with the problem of the preparation of the 
remaining collective mode C-2k in a desired squeezed vac- 
uum state. This is done in what we call the second step of the 
preparation, in which we first adjust the driving lasers to prop- 
agate along the anti-clockwise mode a_ . We then send series 
of pulses of frequencies, phases and amplitudes the same as in 
the above first stage. As a result of the coupling to the cavity 
mode a_, the interaction is now governed by the Hamiltonian 
(l20l . and therefore after the unitary squeezing transformation 
the Hamiltonian of the system takes the form 

He = So(-£o)S±fc(-£i)ff e So(£o)S±fc(£i) 



state >S±fc(£i)|0c 2fc , 0c_ 2fc ) associated with the superposition 
of two, position dependent counter-propagating modes. It is 
interesting to note that only a single set of the laser parame- 
ters, the Rabi frequencies and phases, was required to create 
the pure squeezed states. 

The pure squeezed states are of particular importance in 
quantum physics as the density operator of the system pre- 
pared in the pure squeezed state cannot be factorized into the 
product of density operators of the individual modes. Thus, 
we may conclude that the atomic ensemble, after the interac- 
tion with the sequences of the laser pulses considered in this 
section is prepared in single and two-mode entangled states. 



IV. CREATION OF A FOUR-MODE SQUEEZED STATE 
WITH TWO ATOMIC ENSEMBLES 



It is particularly interesting to generate entangled states be- 
tween two (n = 2) atomic ensembles. Since the cavity is 
composed of two counter-propagating modes, the interaction 
between the collective modes and the pair of cavity modes can 
produce squeezing of the multi-mode variety, with a possibil- 
ity to create a four- mode squeezed state fl8l l48l l49ll . We will 
illustrate how it is possible to create a four-mode squeezed 
state between two atomic ensembles by a separate address- 
ing of the collective modes with sequences of laser pulses of 
suitably chosen Rabi frequencies and phases. 

Let us demonstrate a procedure which constructs the four 
in) 

collective bosonic modes C ±2k in the following pure four- 



mode squeezed state 



V^Si - Pli (a-Cok + a\C-2k + H.c.) . (31) |*) = S ±k (0\0 c w,0 c 



(34) 



As above in the case of the coupling to the cavity mode a + , the 
Hamiltonian (IBTt describes a system of two independent lin- 
ear mixers. Hence, the state of the system will evolve during 
the interaction towards its stationary value, and after s suitably 
long time, ~ 2/k, the transformed system will be found in the 
vacuum state. 

Thus, after the second step of the preparation, the trans- 
formed system is found in the pure vacuum state determined 
by the density matrix of the form 



where the unitary operator has the form 



exp 



°2fc °-2fc 



°-2fc ( - / 2fe 



(2)t r ,(2)t 



C^C. 



-2k 



H. 



•)}• 



(35) 



p( T ~ 4//s) = |*X*I> 



(32) 



where 



I*} = S (-£o)S ±fe (-£i)|*) 

= |0 a+ ,0 a _,0 Cofc ,0c 2fc ,0c_ 2 J (33) 

represents the vacuum state of the transformed system and the 
ket | "J) represents the final stationary state of the system. 

If we now perform the inverse transformation from p to p 
or from |\&) to \ty), we find that in the steady-state the cavity 
modes are left in the vacuum state and the atomic ensemble 
is prepared in single and two-mode squeezed states. In other 
words, we find that the collective mode Cok is prepared in the 
one-mode squeezed vacuum state £o(£o)|Oc k)> whereas the 
collective modes C±2k are in the two-mode squeezed vacuum 



and £ is the squeezing parameter which we take to be a real 
number. As in the above cases of one and two-mode squeez- 
ing, the squeezing parameter can be adjusted and controlled 
by the driving laser fields. 

The action of the operator ( 1351 ) on the vacuum state is to 
produce the four-mode squeezed state with correlations exist- 
ing between counter-propagating field modes of each and dif- 
ferent ensembles. Actually, the state |\&), after a suitable local 
transformation can be related to a weighted four-mode square- 
type cluster state, as shown by Menicucci et al. j49ll . It has 
been demonstrated that the cluster state (f34-b can be realized by 
engineering nonlinear type optical interactions in quasi-phase- 



matched materials 
(PPKTP) crystals 



15011 . such as periodically poled KTiOP04 
1 5111 . As we shall see, the state could be 
experimentally realized with two ensembles of cold atoms 
placed inside a ring cavity. 

Before introducing the procedure, we would like to point 
out that a transformation of the density matrix of the sys- 
tem with the unitary operator (f35T > does not lead to a den- 
sity matrix describing a simple system of independent linear 
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mixers. Therefore, we first transform the state to a new 
state |$) = r|^), with a unitary operator T chosen such 

(n) 

that it transforms the field operators C^k to new operators 
d^m = TC^ k T\ which are linear combinations involving 
the operators of different ensembles only, i.e. 



a 



AC 



VT+X 2 



4 2) 



AC 



(i) 

2k 



C 



(2) 
2k 



Vl + A^ 



AC 



(2) 
-2k 



c 



(1) 

-2k 



,(2) 



c 



(2) 
-2k 



AC 



(i) 

-2k 



Vl + A^ 



Vl + A^ 



(36) 



where A = (l + \/5)/2. 

With this carefully chosen transformation, we find that the 
state |<1>) takes a form 



|$) = exp 
x|f> 



-i ( xc^c^ 



1 



21: 



— -C 



(2)t^(l)t 



H.c. 



.0. 



0, 



. (I 



SI^AO^C-C/^IO^, 



0, 



.0 



$37) 



where S±l(X(,) and S±l(— £/A) are squeezing operators in- 
volving field modes of the ensemble 1 and 2, respectively. As 
we shall see, the advantage of working with the transformed 
state |$) rather than the state |\I>) is that the density matrix 

transformed with the operators S±fc(A£) and S±l(— £/A) de- 
scribes a simple system of a linear mixer of a cavity mode 
with one of the field modes. 

We now proceed to perform the construction of the four- 
mode squeezed state (l34l . which is done in four steps, each 
separately addressing one of the collective modes. Since the 
collective modes are orthogonal to each other, an arbitrary 
transformation performed on one of the modes will not affect 
the remaining modes. 

In the first step, we adjust the driving lasers to propagate 
in the direction of the clockwise mode a + . In this case, the 
dynamics of the system are described by the Hamiltonian ( TT~8b . 
and after a unitary transformation of the density matrix 

S^/X)S^(-\OTpT^S ( il(XOS { il(-C/\) = pi, (38) 



with the squeezing parameter 



^ 2A 



A.1 + Ps; 



and with the laser Rabi frequencies and phases such that 

<A«„ = 4>» n = °) n = l,2, 



(39) 



(40) 



we find that the master equation of the transformed density 
matrix is of the form 



dt 



Pi = -i[H^pi] + C c pi, 



where 



m = 



Ail ( C Ok + ^ C 0k ) + A>2 ( ^ C 0k ' + Q 



■,(2) 



-r(l)t 



.(2)t 
J 0k 



H.c. 



(42) 



We see that the choice of the laser parameters (l40l > results in 
the collective mode Cil , out of the four modes involved in 
the squeezing operator ( T35l l. being effectively coupled to the 
cavity modes, the other modes are decoupled. It may appear 
surprising that only one of the four collective modes is effec- 
tively coupled to the cavity mode. However, this is merely a 
consequence of the form of the effective interaction Hamilto- 
nian < fT~8b . It has an obvious advantage that we can separately 
address each of the collective modes. One may also notice 

from Eq. d42l that apart form the collective mode C^ k \ the 
(n) 

modes Cg fc are also coupled to one of the cavity modes. Since 

(n) 

the modes Cg fc are not involved in the multi-mode squeezed 
state ( 1341 ) whose construction we are interested in, we do not 
consider their evolution. 

Fortunately, the master equation (l40b is of a similar form 
as Eq. (f26b . Therefore, we can follow the same arguments as 

before to conclude that the mode C^jJ will evolve in time and 
definitely after a sufficiently long time, ~ 2/k, the mode will 
be found in a stationary vacuum state. The other three modes 
will remain in undetermined states. 

We now turn off the lasers propagating in the direction of 
the clockwise mode, and perform the second step, in which 
we first adjust the driving lasers to propagate in the direction 
of the anti-clockwise mode a_. In this case, we adjust the 
squeezing parameters such that 



A*i Xf3 U2 , {3 S2 — Xf3 Sl , 
<t>un = 4>a n =0, n = 1,2, 



(43) 



to obtain the same value for the squeezing parameter now 
given by 



2A 



In 



Pu 2 + /3 S1 

P« 2 Psi 



(44) 



For this choice of the parameters, and bearing in mind that 
the dynamics of the system are now governed by the Hamilto- 
nian ( f20b . we find that after a unitary transformation, Eq. d38l ). 
the master equation of the transformed density matrix takes 
the form 



— p 2 = -i[H 2 ,p2] + £cP2, 



(45) 



with 

H 2 = 



Pu 2 ( Co*? -+ 



+ & I (AC$ t + C ( 



.(2)t 
'Ok 



+ V / ( 1 + A 2 ) ~ /%) C { Xa{ + H.c. 



(46) 



We see that the choice of the laser parameters, Eq. d43l , results 

(41) 

in the collective mode C^l k to be effectively coupled to the 
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cavity modes, with the other modes decoupled. The dynam- 
ics of the mode are determined by the master equation d45b . 
which is of the similar form as Eq. ( f40b - Thus, we can use 
the same arguments as in the previous step, and conclude that 



,(2) 



will 



after a sufficiently long time, ~ 2/k, also the mode C_2 k 
be found in a stationary vacuum state. 

In the next step, we change the direction of propagation 
of the driving lasers back to the direction of the clockwise 
mode a_|_. Similarly as in the above two steps, we first perform 
the unitary transformation of the density matrix, Eq. d38l , with 
the laser parameters chosen such that 



p«i — A/3 U2 , (3 S2 — \p si , 

4>ui = 0s i = o, 4> U2 = 4> S2 



and 



2 



P«2 Psi 
HU2 Psi 



(47) 



(48) 



With this choice of the laser parameters (|47| |, we find that now 

(2) 

the collective mode C^J is the only mode coupled to the cav- 
ity modes and the dynamics of the mode follows the same 
pattern as the modes considered in the above two steps. Thus, 
we may conclude that after a sufficiently long time, also the 

(2) 

collective mode C^J will be found in a stationary vacuum 
state. 

Finally, in the fourth step, we prepare the remaining 

mode C_lk m me vacuum state. To do this, we again change 
the direction of the driving lasers to propagate in the direction 
of the anti-clockwise mode and choose the laser parameters 
such that 



and 



Pu2 ^Pui J Psi ^Ps2 1 



2 n \p ui -/3 S2 



(49) 



(50) 



Following the same procedure as in the previous three steps, 
one can easily show that the choice of the laser parameters, 
Eq. j49l , results in a master equation for the transformed den- 
sity operator determined by the interaction Hamiltonian in- 
volving the collective mode C_l k . Thus, as a result of the 
damping of the cavity modes, the collective mode will evolve 
towards a vacuum state. 

In this way, the final state of the transformed system is a 
four-mode vacuum state determined by the density operator 



represents the vacuum state of the transformed system and the 
ket | "J) represents the final stationary state of the system. 

We see that in this practical scheme we can generate on de- 
mand an arbitrary multi-mode squeezed state. Generalizations 
of the scheme to the case of three or more atomic ensembles 
located inside a ring cavity are also possible and are trivial. 



V. CONCLUSIONS 



We have proposed a procedure which constructs multi- 
mode squeezed (entangled) states in ensembles of cold atoms 
located inside a ring cavity. The procedure is referred to 
practical ring cavities composed of three mirrors which cre- 
ate two counter-propagating modes of degenerate frequencies. 
In contrast to the case of hot atomic ensembles, where the 
one-mode approximation can be used to model the dynam- 
ics of the atoms, the cold atomic ensembles are equally cou- 
pled to both of the cavity modes. Using the the Holstein- 
Primakoff representation of angular momentum operators, we 
have expressed atomic systems in terms of orthogonal collec- 
tive bosonic modes and then have demonstrated how to con- 
struct multi-mode squeezed states by a separate addressing 
of the collective modes. In particular, we have shown how 
to construct one and two-mode squeezed states with a sin- 
gle ensemble and four-mode squeezed states with two ensem- 
bles of cold atoms. The procedure constructs squeezed states 
from the vacuum by a unitary transformation associated with 
the collective dynamics of the atomic ensembles subjected to 
driving lasers of a suitably adjusted amplitudes and phases. 
The lasers prepare, with the help of cavity dissipation, the col- 
lective modes in a desired stationary squeezed state. 

The procedures which construct squeezed (entangled) 
states by using only simple unitary transformations, the theme 
Krzysztof Wodkiewicz was interested for many years, offer 
further interesting prospects for the study of multi-mode en- 
tanglement and creation of entangled cluster states. 
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